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Abstract. We consider an auxiliary fields formulation for the general fourth-order 
gravity on an arbitrary curved background. The case of a Ricci-flat background 
is elaborated in full details and it is shown that there is an equivalence with the 
standard metric formulation. At the same time, using auxiliary fields helps to make 
perturbations to look simpler and the results more clear. As an application we 
reconsider the linear perturbations for the classical Schwarzschild solution. We also 
briefly discuss the relation to the effect of massive unphysical ghosts in the theory. 
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1 Introduction 

The long-standing interest to higher derivative gravity theories is based on their critical impor¬ 
tance for quantum extensions of gravity. In this case the general fourth order gravity represents 
a minimal UV completion in the vacuum (purely metric) sector of the theory. Such a theory 
is renormalizable at both semiclassical [1] (see also [2, 3] for the introduction) and quantum 
[4] levels. At the same time the higher derivative quantum gravity can not be considered as a 
completely consistent theory, because its spectrum of states includes unphysical massive ghost 


- the spin-2 state with negative kinetic energy. Starting from 70-ies there were several different 
proposals on how to solve the problem of massive ghosts [4, 5, 6, 7], mainly based on the quan¬ 
tum field theory methods, but all of them did not lead to conclusive results [8]. Very recently 
one of us proposed a more simple classical way of dealing with ghosts [9]. The new approach to 
the problem is based on the study of stability of the classical solutions in the low-energy sector 
of the theory, and in this respect it is similar to the consideration based on the effective quantum 
field theory ideas [10]. 

The linear stability of the given solution of some differential equation is usually sufficient for 
the stability at any perturbative level. Therefore, in order to observe the effect of the spin-2 
ghosts one can simply look at the asymptotic behaviour of the gravitational waves on the given 
background. In the case of cosmological solutions discussed in [9], we could observe that the 
ghost does not produce instability below the threshold, which approximately corresponds to the 
Planck scale of frequencies. Only starting from the Planck scale one can observe an explosive- 
type dynamics of the gravitational waves. The situation is completely different for tachyons, in 
this case an instability takes place independent on the energy scale. This important difference 
was recently discussed in [11] and [12]. 

The stability of the cosmological background due to the huge mass of the ghost [9] is in 
accordance with the previously known result for the particular deSitter background [13, 14, 15]. 
In these papers the anomaly-induced semiclassical contributions were taken into account, but 
finally were shown to be irrelevant [16]. It would be definitely interesting to consider other 
physically relevant background solutions, in particular for the Schwarzschild metric. During 
some time the unique study of the stability of the d = 4 classical Schwarzschild solution in the 
presence of higher derivative gravitational terms was the work by Whitt [17]. The main result 
is that the solution is stable in the case of fourth derivative gravity. Recently there was another 
paper by Myung, [18], which argued in favor of the opposite result, by using the analogy with 
the well-known Gregory-Laflamme instabilities [19] (see also [20]). The difference between [17] 
and [18] does not concern the equations for perturbations, but only their solutions. 

In this short report we discuss the procedure of separating Ricci-flat and higher-derivative 
perturbations of [17]. For this end we employ the method of auxiliary fields, as it was described 
in a recent paper [21] (see also [22] for standard reviews and further references). In our opinion, 
the possibility to present the general fourth derivative gravity in the second-order form by means 
of tensor auxiliary fields has some independent interest. On the top of that we apply this new 
description to the analysis of metric perturbations on the background of Schwarzschild solution. 

2 Standard approach to fourth-derivative perturbations 

In this section we summarize the Whitt approach [17], which will be verified in the next section 
by using auxiliary fields. According to [17], it is possible to explore the stability by means of 
the second-order equations for the perturbation of the Ricci tensor, instead of the fourth order 
equations in the metric perturbations. This approach helps to simplify the standard treatment 
of the stability problem in the theory of the fourth order in derivatives. 
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Let us consider the gravitational action of the fourth order theory, 


s[g,u] = (1) 

where is the Weyl tensor, = E + — 1/2, R^). In this expression we use 

notations which are typical in quantum gravity, in particular the transition to the notations of 
Whitt [17] can be performed as follows: 
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Furthermore, we disregard the Gauss-Bonnet term with the integrand E = ~ 

since it is not relevant at the classical level. 

The variation of the action (1) with respect to the metric yields 
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Taking trace we arrive at 


^ Di? = 0. 


( 4 ) 


According to this equation, the scalar curvature obeys the covariant minimal Klein-Gordon 
equation. Furthermore, in case of a; = 0 the last formula boils down to the non-dynamical 
constraint on the scalar curvature R = D. The choice a; = 0 corresponds to the case of a pure 
Weyl-squared higher derivative part in Eq. (1). The higher derivative part of the action possess 
local conformal invariance, the symmetry which is violated only by the Einstein-Hilbert term. 

Equation (3) shows that any vacuum solution of Einstein equations is also a solution for the 
theory of the fourth order (1). For the perturbations this may be not true, but in the case when 
the background metric satisfies vacuum Einstein equations, they can be classified into the two 
distinct classes: 


• The perturbations which satisfy Einstein’s equations in vacuum, R^uigap + ^o/?) = 0. 

• The ones which do not satisfy these equations. 

In the first case we have exactly the same problem as in Einstein’s theory because 

^Hui9a/i T ^a/s) R/ii/igaf^) T Oj (5) 

where is the first order expansion of the Ricci tensor, which will be specified below. 

This ensures, for example, the stability of the Schwarzschild metric, as it was discussed in 
[23, 24, 25]. In the second case one can consider the perturbation of the Ricci tensor, 

R'hp = Rp.i' + ■ ( 6 ) 
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Substituting (6) into (3) and remembering that = 0, we obtain the equation for the 

linear perturbations of the Ricci tensor 
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( 7 ) 


where we use notation = g^j_yR^jl}. Let us note that in the Ricci-flat case R^^^ is also the 
first order term in the expansion of the curvature of the second type. 

Taking the trace of (7) we get 


ArIL + ^ DRlb =0, 

A 

which is the linearized form of Eq. (4). For a; = 0 this gives the constraint R^^^ = 0. 

In what follows we shall limit our analysis to this particular case^. So eq. (7) becomes 


( 8 ) 


□R« + = 0 • 


(9) 


Eq. (9) determines the dynamics of small gravitational perturbations. It turns out that the 
stability of the background Schwarzschild metrics in the fourth order theory in d = 4 is defined 
by the same equations as in the black string case in d = 5 space-time dimensions [19] and as 
in the bi-metric theory of gravity [20]. As it was noted in [18], in these cases one meets a 
well-known Gregory-Laflamme instability [19]. 


3 Using auxiliary fields 

Let us consider the action with auxiliary fields 

S 2 = j -^R + + I , ( 10 ) 

where is an invertible c-number operator depending only on the metric. The variation of 

this action with respect to the field gives 

4>a/3 = - ^Aap^^uR^'' ■ ( 11 ) 

After using the last relation, the action S 2 on-shell becomes 

^2 = y - ^7? - A • (12) 

Comparing (1) with (12) we obtain 

P, d^u,al3 'y9^iu9a(3 i ( 18 ) 

^For different values of w, the trace free part and the trace part of do not decouple leading to a fourth 
order equation for Fpp}, where Rpt} — — jgp,i,R^^P 
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where 7 = (1 — ti;)/3. Here we use the standard DeWitt notation 

^^u,aj3 '^{,3^laQuf} S^iPSua) ■ 

The operator defined by (13) does not have an inverse if 7 = 1/4, that is for oj 
simply assume^ 7 7 ^ 1/4. Then the inverse operator can be easily obtained with 


— 1 Aa/S^rX 
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and one can check that 
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where 


1 — CJ 

1 — 4a; 


(17) 


In this way we arrive at the explicit form of Eq. (10), 

^2 = 1 d^xV^[ -^R + ^ ■ ( 18 ) 

It is natural to check whether the auxiliary field (pajs satisfied some constraints. First of all, 
without losing generality we can assume that (j)a/3 is a symmetric space-time tensor. Furthermore, 
one can use Eqs. (II) and (13) to derive the Bianchi identity for this field in the form 

where (j) = 4A''g^^. (19) 

There is a complicated question on whether the auxiliary field should satisfy the constraint (19) 
at the quantum level, but since in the present work our attention is restricted to the purely 
classical case, this relation should be respected. 

The action (18) describes a dynamical theory of the two symmetric tensor fields g^i, and 
It is important to note that the auxiliary field is dynamical, because variation with respect 
to the metric produce second derivatives of in the equations of motion. The action (18) is 
second-derivative, but it is dynamically equivalent to the original fourth-derivative action ( 1 ). 
This equivalence means one can map any solution of one theory to some solution of another one 
and vice versa. Indeed, this is true for both background and perturbed solutions, so the action 
(18) can be useful to explore the perturbations. 

In order to illustrate how it works, let us expand this action up to the second order 

52 = SP + ( 20 ) 


in the perturbations of both metric and auxiliary field. The expansions of the fields are defined 
as 


dfJ.U ^ dfJ.U 

II 

+ 

(21) 


II 

+ 

(22) 


^In the special case 7 = 1/4 the formulation via auxiliary tensor field is possible, but this field must be traceless. 
It is not difficult to elaborate this version, but there is no much practical interest to do it. 
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In what follows we shall also use notations for the traces h = and = '0^- One 

can easily check that the first terms in the expansions and satisfy the same Bianchi 
identities as the background quantities, 

= I = (^ - ^) • (23) 

Other relevant quantities are expanded to the second order as follows: 

= \/^(l + 2^ ~ + • • • ) • 

Furthermore, the first-order expansion for the Ricci tensor is 

= -^nV-^V^V./i + ivAVX + ^VAV./r^ (24) 

= 2 — V^Vi//i -|- R^h\u + R^hx^ — 2Ra^j3i,h°‘^'^ 

and the second-order term is 

= ]jh'^^{VrVxh^y + V^V^hrX-Vx'^vhr^-VxyiJihrv) 

+ ^ ^v^/i) {Vxh^u - V^hxiy - Viyhxfj) 

+ + ^Vrhxf.'V'^h^ - ^V^/i^Va/i); . (25) 

The Ricci scalar in the first order is 

R^^^ = -ah- , ( 26 ) 

and in the second-order we meet 

r(2 ) = + V^V,/i - V^Vxhi - VxVX) + V;, 

- ^VxhV^h + ^Va- ^Va• (27) 

Replacing these formulas into the action S 2 (18), we obtain the first-order expansion 

= J + ^ [ _ - 002 )] 

- 2(j)f,aR'ph>^'' - ^ (p^aPuh^"' + ^ Opxpf.uh^"' + - 04>i>) 

- +0^"(VAVX-^V^V,/i-inV)}- (28) 
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Finally, the second-order term is 

+ Y +1 - OH) + 

+ {r'' - 20 ^ h!'^)R^R + + 2<P^ph^^hl - ct>^ph^^h)R^, 

- \ ^ ^ 

+ ^ - ^ ^] } • ( 29 ) 

The expression (29) contains the complete equations for the perturbations. However, for the 
sake of simplicity we assume that the background metric satisfies Einstein equations in vacuum, 
such that we take = 0 and ipfiu = 0. After discarding surface terms, we obtain the simplified 
expression 

+ r''[- \^K. - \y^yuh + V^Vxht - Rf^augh^^ 

The equations that determine the dynamics of the perturbations have the form 

- ^ {d^,u,ap - Ogf.uQapW^ = R^i!} (31) 

and 

^ {gfiu^h - g^vVaViBh"^^ + V^Vxh^ + V^Vxh^ - V^V,,/r - - 2R^avgh°‘'^) 

+ V^VaV^;) -F Vi.VaV’^ - - 5'Mi^VaV^V’"^ - 2R^augi^°‘'^ = 0. (32) 

By means of relations (23) and (26) the last equation can be presented in the form 

Y ~ \ = \ “ \ ■ (33) 

Inverting Eq. (31) and inserting it into Eq. (33) one can reproduce Eq. (7). For the 
conformal fourth-derivative case, when cu = 0, the equation is 

□4)] + 2R^^xuR^^^^^ + = 0 , (34) 

where we used the constraint R^^'> = 0 (consequently if) = 0 too). Furthermore, one can impose 
the TT-gauge on the metric fluctuations, 

h = 0, = 0, (35) 
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so that Eq. (24) which relates metric and Ricci fluctuations becomes 


— Cl) (36) 

or, in terms of the auxiliary field, 

nhf,^ + 2Rrf,xuh^^ = (37) 

In the last formula we took into account the constraint = 0 which holds in the case a; = 0. 
The same constraint, together with the TT-gauge (35), reduce the Eq. (32) to the relation 

ahf,iy + = -K^ • (38) 

Taken together, Eqs. (37) and (38) lead to the massive equation for the auxiliary field alone, 

~ 6 • (39) 

AC 


4 Brief discussion 


In the previous section we have shown that the two descriptions of perturbations are indeed 
equivalent. Let us see how one can apply the auxiliary fields description for the black-hole 
background case. 

Starting from equations (34) and (36) one can define a new field 


K 


}_h — —h _ p(i) 

tt A 


which satisfies the massless equation (on a Ricci-flat background) 


(40) 


= 0 . 


(41) 


This is exactly the equation for the Ricci-flat perturbations in Einstein theory, which were 
explored by Regge, Wheeler et al in [23, 24, 25]. It is known to have no unstable modes. 
Therefore, as far as unstable modes are concerned we can set = 0, i.e. = —h/Miz/K^, 

making Eq. (36) to take the same form as Eq. (34). 

In order to establish the relation between massive ghosts and instabilities, one has to note 
that the ghost with the mass of the Planck order of magnitude does not produce instabilities in 
flat space or in a very weak gravitational field [11, 10]. The practical consequence for the black 
hole case should be a selection rule for the boundary conditions of the perturbations at space 
infinity. Neglecting the oscillatory part, we can impose the requirement for those perturbations 
which are related to the ghost instabilities that 'ipfj.u, h^u vanish for large values of r, where the 
background solution approaches Minkowski space. Only those solutions for the perturbations 
should be permitted, if we intend to estimate the effect of massive ghost on the instabilities on 
the black-hole background. 



The stability of the perturbations is completely defined by eq. (34). This equation is known 
to possess a spherically symmetric unstable mode provided 

0 1 

VP ^ '^s 

where rs is Schwarzschild radius [19]. Such a mode grows exponentially in time as > 0, 

spatially vanishes at infinity and is regular at the future horizon. This mode is present in both 
(as also noted in [18]) and h^,^. 

The Gregory-Laflamme unstable mode exponentially decays in the radial direction, leading 
to 

(V> ~ (42) 

Eq. (42) tells us that in any given finite space point at some instant of time the perturbations 
become large and then show an unrestricted growth. According to our previous considerations, 
one can consider this as an indication that the Gregory Lafiamme instability is not related to the 
presence of massive ghosts in the spectrum. Finally, we note that the end-state of this instability 
is thought to be a black hole with a massive graviton halo 

5 Conclusions 

We presented a new form of the action and of equations of motion for the higher derivative 
gravity, which uses auxiliary fields on an arbitrary curved space-time background. Previously, 
the same representation has been known in flat space [27] and for the special de Sitter background 
[28]. The main result is the eq. (29) from which one can easily extract the full set of equations 
in the physically interesting cases or approximations. For the sake of simplicity, we elaborate 
only the perturbations for the Ricci-flat background and arrive the coupled set of equations for 
metric and auxiliary held perturbations. 

As an application we consider linear perturbations around the static spherically symmetric 
black hole. The use of auxiliary tensor held enables one to distinguish the perturbations related 
to the massive unphysical ghost which is present in the spectrum of the theory. As a result we 
can see that these perturbations can not be identihed with the Gregory-Lahamme instability. 
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